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Abstract 

Starting from a master equation in a quantum Hamilton form we study ana- 
lytically a nonequilibrium system which is coupled locally to two heat bathes 
at different temperatures. Based on a lattice gas description an evolution 
equation for the averaged density in the presence of a temperature gradient is 
derived. Firstly, the case is analysed where a particle is removed from a heat 
bath at a fixed temperature and is traced back to the bath at another temper- 
ature. The stationary solution and the relaxation time is discussed. Secondly, 
a collective hopping process between different heat bathes is studied leading 
to an evolution equation which offers a bilinear coupling between density and 
temperature gradient contrary to the conventional approach. Whereas in case 
of a linear decreasing static temperature field the relaxtion time offers a con- 
tinuous spectrum it results a discrete spectrum for a quadratically decreasing 
temperature profile. 
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I. INTRODUCTION 



Within the classical equilibrium thermodynamics the Carnot cycle had been considered 
as the standard example for studying the efficiency of heat devices acting between two heat 
bathes. The analysis is restricted to reservoirs with fixed temperatures and furthermore 
to the reversible limit. Consequently Carnot 's engine works only with a fixed rate of heat 
provided from the reservoirs. 

Recently, Velasco et al0 studied a finite time Carnot refrigerator to get an upper bound for 
the coefficient of performance of endoreversible refrigerators. They found an upper bound 
for the mentioned coeffcient depending on the ratio between the temperature of the cold 
and the hot reservoir. 

There is a general interest in nonequilibrium systems with two temperature^H^. The analysis 
is motivated by searching for some generic features of nonequilibrium steady states. In par- 
ticular, the question appears for a universal behavior under nonequilibrium conditions. As 
an example, a two-temperature, kinetic Ising model is investigated^ extended to a diffusive 
kinetic system i ni. The authors found a bicritical point where two nonequilibrium critical 
lines meet. The analysis is strongly supported by Monte Carlo simulations in two dimen- 
sions. Recently, a similar simulation has been performed studying a two-temperature lattice 
gas model with repulsive interactions^. A complete different approach is used ini where a 
thermally driven ratchet is studied under periodic, dichotomous temperature changes. The 
behavior of the engine is significant different from a quasistatically working one. InS the 
author used a local heat conduction operator to study the corresponding thermal processes 
observed in complex fluids. Another approach consists of the analysis of a cyclic working 
thermodynamic devise driven by an external applied steady flowi. A nonlinear oscillator 
coupled to various heat batheshad been considered as a simple toy modell. 
Here we are interested in a 'local Carnot engine' on a lattice gas, i.e. each point of a lattice 
will be contacted with two heat bathes at local different temperatures. A particle taken away 
from a reservoir is created at a lattice point i where the creation rate depends on the local 
temperature related to this lattice point. In the same way a particle is annihiliated from the 
neighboring point j. This particle is removed to a reservoir on a different temperature. As 
the result we consider the hopping of particles from a lattice point to its neighboring one 
wheras both points are in contact to heat bathes on different temperatures. Alternatively, 
a model is studied where a particle at a certain lattice point is able to change its state, 
may be from spin up to spin down, however this ffip process is organized by coupling to 
local bathes. As before the up and down states are also coupled to reservoirs at different 
temperatures. As a useful method to study such situations we apply the quantum formalism 
for nonequilibrium processes based upon spin variables. 



II. QUANTUM APPROCH TO NONEQUILIBRIUM 

The analysis is based on a master equation 

dtPin,t)=L'P{n,t) (1) 

where P{n,t) is the probability that a certain configuration characterized by a state vector 
n = (ni, 77-2 . . . ■n.Af) is realized at time t. In a lattice gas description each point is either 
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empty or single occupied nt = 0, 1. These numbers can be considered as the eigenvalues of 
the particle number operator. The dynamics is determined completely by the form of the 
evolution operator L', specified below, and the commutation relations of the Pauli-operators. 
Thus, the problem is to formulate the dynamics in such a way that this restrictions in the 
occupation number are taken into account explictly. The situation in mind can be analyzed in 
a seemingly compact form using the master equation in a quantum Hamilton formalismSHli, 
for a recent review seeElS. Within that approach^ the probability distribution P{n,t) is 
related to a state vector | F{t)) in a Fock-space according to P{n,t) = {n \ F{t)). The basic 
vectors | n) are composed of Pauli-operators. Using the relation 

\F{t)) = Y.P{n,t)\n) (2) 

rii 

the master Eq. (|T]) can be transformed into an equivalent one in a Fock-space 

d, I Fit)) = L I Fit)) (3) 

where the operator L' in ([l|) is mapped onto the operator L in Eq.(|^). It should be empha- 
sized that the procedure is up to now independent on the realization of the basic vectors. 
Originally, the method had been applied for the Bose caseli3llilEl. Recently, an extension 
to restricted occupation numbers (two discrete orientations) was proposedc3il3. Further 
extensions to p-fold occupation numbersii as well as to models with kinetic constraints are 
possible!!. 

As shown by Doi0 the average of an arbitrary physical quantity Bin) can be calculated by 
the average of the corresponding operator Bit) 

(Bit)) = Y.Pin,t)Bin) = {s \ B \ Fit)) (4) 

rii 

with the state function {s \= \- The evolution equation for an operator Bit) reads now 

d,{B) = {s\[Bit),L]\Fit)) (5) 

As the result of the procedure, all the dynamical equations govering the classical problem 
are determined by the structure of the evolution operator L and the commutation rules of 
the operators. In our case the dynamics will be realized either by spin-flip or by exchange 
processes, respectively. 

III. COUPLING TO HEAT BATHES 

The evolution operator for a local flip-process readJ^^ 

Li = A(4 - didl) + 'jidi - d\di) (6) 
where A and 7 are independent flip-rates. 

The occupation number operator ni = d\di is related to the spin due to the relation Si = 
1 - 2n,-. 
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A generalization to processes under the influence of a heat bath with a fixed temperature T 
is discussed inii. As demonstrated inSil the evolution operator has to be replaced by 

L = z/^ [(1 - di) exp{-f3H/2)dlexp{(3H/2)] 

+ [(1 - 4) exp{-pH/2)d,exp{PH/2)] (7) 

where i/ is a the flip-rate defined on a microscopic time scale; f3 = is the inverse tem- 
perature of the heat bath and H is the Hamiltonian responsible for the static interaction. 



A. Flip-dynamics 

Whereas by Eq.(|^) the coupling to a global heat bath is realized we discuss now a further 
generalization by introducing two local heat bathes with different temperatures T and T', 
respectively. The two reservoirs are coupled directly to each lattice point. This situation 
can be described by an evolution operators 



Here /i is a characteristic energy which is necessary to remove a particle from the heat bath 
or to give it back to the bath. The approach reminds of using a grand canonical ensemble 
in equilibrium statistics. Therefore, the quantity /i plays the role of the chemical potential 
assumed to be identically for both processes under consideration. The potential fi can be 
positive or negative. Taking into account that the occupation operator rii has the eigenvalues 
or 1 we get 

Thus, the operator di annihiliates a particles at the temperature T independently on the 
temperature T' of the other bath. Contrary, the operator dl creates a particles at the 
temperature T'. The evolution operator describes the process of annihiliation and 
creation of particles within the system at different temperatures. Using the Eq.(|^) and 
the algebraic properties of Pauli-operators, the evolution equation for the averaged densitiy 
reads 

iy-'dt{ni) = exp(-^/(2r'))(l - n,) - exp(^/(2r))(n,) (10) 
This equation can be solved easily. It results a stationary state at an effective temperature 



(n)s = , , with — = -(- + — ) (11) 

^ ' l + exp(/i/re) Te 2^T ^ ' 



In a spin representation we obtain 



eM/2r _ g-M/2r' 
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In the special case that T = T' the stationary solution coincides with the conventional 
equilibrium solution 

{S)s=i^n\.^ {n)s = ^;:j^^ (13) 

If the temperature of one of the heat bathes tends to infinity (for instance T' —>■ oo) the 
stationary solution is 

{S)s = tanh^ 

When both temperatures T and T' are infinitesimal different from each other T' = T + AT 
the averaged occupation number is 

The Fermi-distribution as the equilibrium solution is modified in lowest order in AT by an 

additonal term proportional to the Bose-distribution. 

The relaxation time r related to the Eq.(|T^) is simply given by 

(^^)-^ = exp(^)+exp(-^) (15) 

The relaxation time for T' ^ T is either enhanced for T' < T or diminished in the opposite 
case. 



B. Exchange Process 

Up to now we have analysed the case of independent flip processes (annihiliation-creation- 
processes) at different temperatures without an internal coupling between the active parti- 
cles. In the following, we discuss the situation that the particles can exchange their mutual 
position; with other words hopping processes are allowed between neigbored sites under the 
influence of the coupling to local heat bathes. The evolution operator reads 



<ij> 



did\)e ^^^^^^U\d,e^^^^ 



(16) 



It describes the exchange process between two adjacent neighboring sites, where the lattice 
site i is coupled to the bath at the temperature Tj and the site j is related to T,-, repectively. 
The evolution equation for the averaged density can be written in the form 

y-'dt{nr) = Y(^j)exp[-i^/2)i^ - h] - (n,) exp[-(/i/2)(l - h] 

j(r) ^ i i ^ 

-2(n,n,)sinh[(/x/2)(l-l)] (17) 

ij ij. 

In the special case of fixed temperatures Tj = Tr = T the last equation is reduced to the 
conventional diffusion equation in a discrete representation. Here, the case will be studied 
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assuming a small temperature gradient. Moreover, Eq.(|T^ is investigated in the continuous 
limit leading to the evolution equation for the density n{x, t) = {n^{t))l~'^ (we set the lattice 
size / = 1) 

y-^dtn = y^n + iiny^^ + li{yn-V{^)^ (18) 



To derive this equation we have neglected the bilinear terms in Eq.(|T^ which give only rise 
to higher order corrections in the density. Due to the conservation of the spins within the 
exchange model the evolution equation can be rewritten as a continuous equation with the 
current j{x, t) 

J = -uVn - ui2nV{l/T) (19) 

In contrast to the conventional nonequilibrium thermodynamics a bilinear coupling between 
the density and the temperature gradient is included in the current. Such a nonlinear 
coupling may change the physical behavior. 
Using natural boundary conditions the stationary solution is 

n(f) =noexp(--^) (20) 
T{x) 

That means, the local density is determined by the local temperature in accordance with 
the local temperature attached to each lattice site. 
Let us consider the special case assuming that 

V(^) = 2c 

where c is a constant vector leading to a decreasing temperature profile 

T{x) = J\ ^ ^ (21) 
2To(c ■ x) + 1 

where Tq is an arbitrary initial temperature. Eg . (pISD can be solved making the ansatz 
n{x,t) = ^{x)ip{x,t). Chosing 

$ = $oGxp[— |u(c • x)] 

ipix, t) obeys 

^ = -[-zyV^ + z/(^c)V (22) 

which is nothing else as the Schrodinger equation for a free particle. The relaxation time is 
the inverse eigenvalue of the Hamiltonian H = — i/V^ + z^(/ic)^ defined in Eq. (p2D : 

_ 1 _ 1 

~ i/(P + (/xc)2) ~ z/(P + /i2(V(l/T))74) ^^^^ 

where k is the wave vector. There is a gap in the quasi-continuous relaxation spectrum for 
k 0. Moreover, the relaxation time depends on the temperature gradient. The solution 
for the density is 
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n{x, t) = uq exp 
Let us consider the further solvable case 



^—jJLC + IK) ■ X 



(24) 



V(^) = 2bx 

which leads to a quadratically decreasing temperature profile, b is a constant. 

To 



T x) 



Tobx^ + 1 



The same procedure as used before yields the Hamiltonian H of the d-dimensional harmonic 
oszillator 

H = -z/V^ + v{^ihfx^ - ufibd 

It results a discrete relaxation time spectrum where the ground state energy of the harmonic 
oscillator is cancelled due to the last term in H. 

r = TT ^ (25) 

2z//i6(mi + 1712 + ... + rrid) 

where the are integer numbers. Obviously, the analysis can be extended to other static 
temperature profiles such as an arbitrary radial symmetric one without changing the results 
substantially. 

A temperature flow is allowed if the heat bathes are coupled. This process leads to an 
equalization of temperatures. Assuming that this process follows the conventional heat 
conduction equation 

1 f2 
TixA) = - — — — 77-exp( — ) A heat conduction 

we conclude that Eq.(0) can also be transformed into a Schrodinger-like equation 

dttl){x, t) = - + V{x, t)] t) (26) 

where the potential V{x, t) is known. The procedure yields a similar result for an arbitrary 
temperature distribution of the form T{x,t) = t~°'g{x'^ /t). 



IV. CONCLUSIONS 

Here we have considered a generalization of the well known Carnot cycle with open flow. 
Each point of a lattice is related to a local heat bath hold on different temperatures. Intro- 
ducing the conventional temperature means that the system is not too far from equilibrium. 
Consistently with this assumption is the consideration of small gradients in temperature 
leading to a heat transport. Such a temperature gradient is coupled to the creation and 
annihiliation of particles or to an exchange process where a particle is created at a certain 
lattice point at a fixed temperature and annihiliated at another point with another but 
fixed temperature. Using a quantum formalism for the master equation we can derive an 
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evolution equation for the density which offers aheady in a mean field like approximation 
a bilinear coupling between density and temperature gradients. This leads to a stationary 
state where the local density is related to a local temperature. The dynamics is studied for 
some special cases with a static temperature profile. The relaxation time spectrum offers 
a different behavior depending on the realization of the temperture field. Moreover, the 
relation time depnds also on the temperature field. 
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